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Q: 

This paper participates in the development of a unified approach to the control of aerial vehicles with extended flight envelopes. 

J~~i" More precisely, modeling for control purposes of a class of thrust-propelled aerial vehicles subjected to lift and drag aerodynamic 

\ forces is addressed assuming a rotational symmetry of the vehicle's shape about the thrust force axis. A condition upon 

•^fy , aerodynamic characteristics that allows one to recast the control problem into the simpler case of a spherical vehicle is pointed 

• ■ out. Beside showing how to adapt nonlinear controllers developed for this latter case, the paper extends a previous work by 

Y\ the authors in two directions. First, the 3D case is addressed whereas only motions in a single vertical plane was considered. 

i i Secondly, the family of models of aerodynamic forces for which the aforementioned transformation holds is enlarged. 
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H ' 1 Introduction 

oj ■ 

■ Feedback control of aerial vehicles in order to achieve 
\ some degree of autonomy remains an active research do- 

■ main after decades of studies on the subject. The com- 
!* ' plexity of aerodynamic effects and the diversity of fly- 

. £^ \ ing vehicles partly account for this continued interest. 

■ Lately, the emergence of small vehicles for robotic appli- 
' cations (helicopters, quad-rotors, etc) has also renewed 

5^ . the interest of the control community for these systems. 
Most aerial vehicles belong either to the class of fixed- 
wing vehicles, or to that of rotary- wing vehicles. The first 
class is mainly composed of airplanes. In this case, weight 
is compensated for by lift forces acting essentially on the 
wings, and propulsion is used to counteract drag forces 
associated with large air velocities. The second class con- 
tains several types of systems, like helicopters, ducted 
fans, quad-rotors, etc. In this case, lift forces are usu- 
ally not preponderant and the thrust force, produced by 
one or several propellers, has also to compensate for the 
vehicle's weight. These vehicles are usually referred to 
as Vertical Take- Off and Landing vehicles (VTOLs) be- 
cause they can perform stationary flight (hovering). On 
the other hand, energy consumption is high due to small 
lift-to-drag ratios. By contrast, airplanes cannot (usu- 
ally) perform stationary flight, but they are much more 
efficient energetically than VTOLs in cruising mode. 



Control design techniques for airplanes and VTOLs have 
developed along different directions and suffer from spe- 
cific limitations. Feedback control of airplanes explicitly 
takes into account lift forces via linearized models at low 
angles of attack. Based on these models, stabilization is 
usually achieved through linear control techniques [15]. 
As a consequence, the obtained stability is local and dif- 
ficult to quantify. Linear techniques are used for hover- 
ing VTOLs too, but several nonlinear feedback methods 
have also been proposed in the last decade to enlarge the 
provable domain of stability [3] [5] [7] [9] . These meth- 
ods, however, are based on simplified dynamic models 
that neglect aerodynamic forces. For this reason, they 
are not best suited to the control of aerial vehicles mov- 
ing fast or subjected to strong wind variations. Another 
drawback of the independent development of control 
methods for airplanes and VTOLs is the lack of tools for 
flying vehicles that belong to both classes. These vehi- 
cles are usually referred to as convertible because they 
can perform stationary flight and also benefit from lift 
properties at high airspeed via optimized aerodynamic 
profiles. The renewed interest in such vehicles and their 
control reflects in the growing number of studies devoted 
to them in recent years [2] [4] [10] [11], even though the 
literature in this domain is not much developed yet. One 
of the motivations for elaborating more versatile control 
solutions is that the automatic monitoring of the deli- 
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cate transitions between stationary flight and cruising 
modes, in relation to the strong variations of drag and 
lift forces during these transitions, remains a challenge 
to these days. A first step in this direction consists in 
taking into account drag forces that do not depend on 
the vehicle's orientation [6], as in the case of spherical 
bodies. 

The present paper essentially aims at extending [6] by 
taking lift forces into account and extending to the 3D 
case a previous contribution [13] concerning vehicles 
moving in the vertical plane (2D case) which shows how, 
for a particular class of models of lift and drag aerody- 
namic forces acting on a wing, it is possible to bring the 
control problem back to the simpler one of controlling a 
spherical body subjected to an orientation-independent 
drag component solely. One can then apply the non- 
linear control schemes proposed in [6] for which quasi 
global stability and convergence results are established. 
The results here reported thus constitute a contribution 
to setting the principles of a general nonlinear control 
framework that applies to many aerial vehicles evolv- 
ing in a large range of operational and environmental 
conditions. 

The paper is organized as follows. In the background 
Section 2, after specifying the notation used in the paper 
and the actuation framework considered to control the 
vehicle, general dynamics equations are recalled and mo- 
tivations for the present study are discussed in relation 
to the limitations of previous results addressing the case 
of spherical bodies. Original results concerning the mod- 
eling of aerodynamic forces acting on symmetric bodies 
and the characterization of a family of models that allow 
one to recast the control problem into the simpler case of 
a vehicle subjected to only an orientation-independent 
drag force are reported in Section 3. Members of this 
family are singled out and tuned by using experimen- 
tal data borrowed from [8] and [14] for elliptic-shaped 
and missile-like bodies. To illustrate the usefulness of 
these results at the control design level with an example, 
Section 4 gives the adapted version of a velocity con- 
trol scheme proposed in [6]. The concluding Section 5 
offers complementary remarks and points out research 
perspectives. 

2 Background 

2.1 Notation 

• The i t h component of a vector x is denoted as xi . 

• For the sake of conciseness, (xii-\-X2j-\-xsk) is written 
as k)x. 

• S{-) is the skew-symmetric matrix- valued operator as- 
sociated with the cross product in R 3 , i.e. such that 
S(x)y = xxy, V(x, y) e R 3 x R 3 . 

• {ei, e2, es} is the canonical basis in R 3 . 



• m is the mass of the vehicle, assumed to be constant, 
and G m is the body's center of mass. 

• X = {O; To, jo ,/co} is a fixed inertial frame with re- 
spect to (w.r.t.) which the vehicle's absolute pose is 

measured, and B = {G; T, J, k} is a frame attached to 
the body. Observe that G and G m may not coincide. 

• The body's linear velocity is denoted by v = 
jf t OG m = (T , jo,k )x = (i,j,k)v. 

• The linear acceleration vector is a = ^v. 

• The body's angular velocity is Co = k)uj. 

• The vehicle's orientation w.r.t. the inertial frame is 
represented by the rotation matrix R. The column 
vectors of R are the vectors of coordinates of i, j, k 
expressed in the basis of 1. 

• The wind's velocity vector v w is assumed to be the 
same at all points in a domain surrounding the vehicle, 
and its components are defined by v w = (i,j,k)v w . 

The airvelocity v a =(h J-> k)v a =(to, Jo> ko)x a is defined 
as the difference between v and v w . Thus, v a = v — v w . 

2.2 Vehicle's actuation 

To cover a large number of actuation possibilities as- 
sociated with underactuated aerial vehicles, and work 
out general principles applicable to many of them, one 
must get free of actuation specificities and concentrate 
on operational common denominators. This leads us to 
assume, as in [6], that the vehicle's means of actuation 
consist of a thrust force T along a body-fixed direction, 
and a torque Yq which allows one to modify the body's 
instantaneous angular velocity Co at will. In practice, this 
torque is produced in various ways, typically with sec- 
ondary propellers (VTOL vehicles) , rudders or flaps (air- 
planes), control moment gyros (spacecrafts), etc. The 
latter assumption implicitly implies that the torque cal- 
culation and the ways of producing this torque can the- 
oretically be decoupled from high-level control objec- 
tives. The corresponding requirement is that "almost" 
any desired angular velocity can physically be obtained 
"almost" instantaneously. Under these assumptions, the 
control of the vehicle relies upon the determination of 
four input variables, namely the thrust intensity and the 
three components of uj. The following complementary 
assumption about the thrust force is made. 

Assumption 1 The thrust force T is parallel to the vec- 
tor k, i.e. T = — Tk withT denoting the thrust intensity. 

The minus sign in front of the equality's right-hand side 
is motivated by a sign convention, also used in [6]. 

2.3 Vehicle's dynamics 

The external forces acting on the body are composed of 
the weight vector mg and the sum of aerodynamic forces 
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denoted by F a . In view of Assumption 1, applying the 
fundamental theorem of mechanics yields the following 
equations of motion: 

ma = rng + F a — Tk, (1) 
= (2 x (ij,k), (2) 

with T and uo the system's control inputs. The motion 
equation (1) points out the role of the aerodynamic force 
F a in obtaining the body's linear acceleration vector a. 
It shows, for instance, that to move with a constant ref- 
erence velocity the controlled thrust vector Tk must be 
equal to the resultant external force F ext := mg + F a . 
When F a does not depend on the vehicle's orientation, as 
in the case of spherical bodies subjected to orientation- 
independent drag forces only, the resultant external force 
does not depend on this orientation either. The control 
strategy then basically consists in aligning the thrust di- 
rection k with the direction of F ext (orientation control 
with uo) and in opposing the thrust magnitude to the 
intensity of F ext (thrust control with T). The almost- 
globally stabilizing controllers proposed in [6] illustrate 
this strategy. However, the production of lift and drag 
forces that depend on the vehicle's orientation may sig- 
nificantly complexify this strategy. In particular, the re- 
sultant force F ext being now orientation-dependent, the 
existence and uniqueness of the equilibrium in terms of 
the vehicle's orientation is no longer systematic, and the 
stabilization of such an equilibrium, when it is locally 
unique, can be very sensitive to thrust orientation vari- 
ations. As a matter of fact, the capacity of calculating 
the direction and intensity of F a at every time-instant - 
already a quite demanding requirement- is not sufficient 
to design a control law capable of performing well in 
(almost) all situations. Knowing how this force changes 
when the vehicle's orientation varies is needed, but is 
still not sufficient. An original outcome of the present 
study is precisely to point out the existence of a generic 
set of aerodynamic models that allow one to recast the 
control problem into the one of controlling a spherical 
body for which strong stability and convergence results 
can be demonstrated. Of course, the underlying assump- 
tions are that these models reflect the physical reality 
sufficiently well and that the corresponding aerodynamic 
forces can be either measured or estimated on line with 
sufficient accuracy. 

3 Modeling of aerodynamic forces 

3. 1 Models of lift and drag forces 

Working out a functional model of aerodynamic forces 
from celebrated Navier— Stokes nonlinear partial differ- 
ential equations governing the interactions between a 
solid body and the surrounding fluid is beyond the au- 
thors domain of expertise, all the more so that spatial 



integration of these equations over the shape of an ob- 
ject does not yield closed- form expressions except in very 
specific cases. Notwithstanding the delicate and com- 
plex issues associated with turbulent flows -a side effect 
of which is the well known stall phenomenon- for which 
no general complete theory exists to our knowledge. We 
thus propose to take here a different route by combin- 
ing a well-accepted general expression of the intensity of 
aerodynamic forces with geometric considerations based 
on the body's symmetry properties. To be more precise, 
let Fd and Fl denote the drag and lift components of 
Fa, i-e- 

F a '= F L + Fd, 

with, by definition, Fl orthogonal to v a and Fd par- 
allel to v a . Consider also a (any) pair of angles (a, j3) 
characterizing the orientation of v a with respect to the 
body frame. The Buckingham ir— theorem [1, p. 34] as- 
serts that the intensity of the static aerodynamic force 
varies like the square of the air speed \v a | multiplied by a 
dimensionless function C(-) depending on the Reynolds 
number R ei the Mach number M, and (a, /?), i.e. 

k a := ^, \F a \ = k a \v a \ 2 C(R e ,M,a,f3), (3) 

with p the free stream air density, and £ an area ger- 
mane to the given body shape. Then, further assuming 
that the direction of F a depends upon the airspeed mag- 
nitude \v a \ via the C(-) function variables (i? e ,M) only 
and that this force does not (or little) depend(s) upon 
the angular velocity cJ, one shows that this theorem in 
turn implies the existence of two dimensionless functions 
Cd(') and Cx(-), and of a unit vector- valued function 
f(-) characterizing the direction of the lift force w.r.t the 
body frame, such that 

F L = k a \v a \C L (R e ,M,a,f3)r(a,f3) x v a , (4a) 
F D = -k a \v a \C D (R e ,M,a,f3)v a , (4b) 
f(a,P).v a =0, (4c) 

In the specialized literature Cd (£ I^ + ) and Cl (g M) are 
called the aerodynamic characteristics of the body, and 
also the drag coefficient and lift coefficient respectively. 

3.2 Aerodynamic models for symmetric bodies 

The expressions (4) of the lift and drag forces hold in- 
dependently of the body's shape, since they are derived 
without any assumption upon this shape. However, 
shape symmetries of aerial vehicles -as well as of marine 
and ground vehicles- are not coincidental. Simplification 
and cost reduction of the manufacturing process, despite 
their importance, are clearly not the main incentives ac- 
counting for the ubiquitous use of symmetric shapes. In 
this respect, Nature was first to give the example with 
most of the animals populating the Earth. On the basis 
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Figure 1. Symmetric and bisymmetric body shapes. 

of this observation, one could figure out numerous prac- 
tical advantages resulting from symmetry properties. 
For flying purposes, however, not all symmetries are 
equally interesting. For instance, the sphere which rep- 
resents the simplest most perfect symmetric 3D-shape is 
not best suited for energy-efficient long-distance flights 
because it does not allow for the creation of lift forces 
which can counteract the effects of gravity, in the same 
way -and almost as well- as wheel-ground contact re- 
action forces for terrestrial vehicles, and buoyancy for 
marine and underwater vehicles. We here consider other 
kinds of symmetries in order to point out aerodynamic 
properties induced by them and their practical interest. 
More specifically, we focus on vehicles whose external 
surface S is characterized by the existence of an or- 
thonormal body frame B = {G; i^j^k}, with k denoting 
the thrust direction according to Assumption 1, that 
satisfies either one of the following assumptions. 

Assumption 2 (Symmetry) Any point P G S trans- 
formed by the rotation of an angle about the axis Gk, 
i. e. by the operator defined by 

0*(O = rot G g(0)(.), 
also belongs to S, i.e. ge(P) G S. 

Assumption 3 (Bisymmetry) Any point P G S 
transformed by the composition of two rotations of angles 
and 7r about the axes Gk and Gj, i.e. by the operator 
by 



ge(-) = (rot G g(0) orot GJ <7r))(-), 
also belongs to S, i.e. ge(P) G S. 



The operator roto$(£)(P) stands for the rotation about 
the axis Ov by the angle £ of the point P. Examples of 
"symmetric" and "bisymmetric" shapes satisfying these 
assumptions are represented in Figure 1. 



For symmetric shapes, i.e. such that Assumption 2 holds 
true, one can define a G [0, tt] as the angle of attac]P~] 

between — /c and v a , and /3 G (— 7r, tt] as the angle between 
the unit frame vector Tand the projection of v a on the 
plane {G; i, j} (see Fig. 2). Observe that this assumption 
also implies that: 

PI : the aerodynamic force F a does not change when 
the body rotates about its axis of symmetry Gk; 
P2 : F a G span{/c,i7 a } . 

Property PI in turn implies that the aerodynamic char- 
acteristics do not depend on /?, whereas Property P2 im- 
plies that i) f is orthogonal to fc, ii) f is independent 
of a, and in) the lift coefficient is equal to zero when 
a = {0,7r}. Subsequently, the expressions (4) of the lift 
and drag forces specialize to 

F L = k a \v a \C L (R e ,M,a)r(f3) x v a , (5a) 
F D = -k a \v a \C D (R e ,M,a)v a , (5b) 
f(f3) = - sin(/3)r+ cos(/3)f. (5c) 

Under the stronger Assumption 3, i.e. when the body's 
shape is also 7r-symmetric w.r.t. the Gj axis, the aero- 
dynamic characteristics Cl and Cd must be tt— periodic 
w.r.t. a. The aforementioned choice of (a, 0) implies that 



cos 



■■p7 ) , P = atan2(v a2 , v ai ), (6) 



and 



r v ai = |v |sin(a)cos(/3), 
i 2 = |v a |sin(a)sin(/3), 
is= ~ \v a \cos(a). 



(7) 



From the definitions of a and r(/3), one then verifies that 



r(/3) x v a = - cot(a)v a r~^k, 

sm(a) 



so that F a = F L + F D becomes 



-kaWa\ 



(C D (-) + C L (.) CO%)Va + ^r\\Va\k 

\ J sin a 



sin (a) 



3.3 Aerodynamic models yielding spherical equivalency 

For constant Reynolds and Mach numbers the aerody- 
namic coefficients depend only on a. By using the above 
relationship, it is a simple matter to establish the fol- 
lowing result. 



1 The angle of attack a so defined does not coincide with 
the one used for airplanes equipped with flat wings which 
break the body's rotational symmetry about Gk [15, p. 53]. 
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Figure 2. Aerodynamic forces and (a, j3) angles. 

Proposition 1 Consider a symmetric thrust-propelled 
vehicle. Assume that the aerodynamic forces are given 
by (5) and that the aerodynamic coefficients satisfy the 
following relationship 

C D (a) + C L (a)cot(a) = C Do , (8) 

with Cd q denoting a constant number. Then the body's 
dynamic equation (1) may also be written as 

ma = rng + F p — T p k, (9) 

with 

T p = T + k a \v a \ 2 ^^-, (10a) 

F p = -k a C Do \v a \v a . (10b) 



with Co and c\ two real numbers, satisfy the condition (8) 
with C D = cq-\-2c\. The equivalent drag force and thrust 
intensity are then given by 

F p = - k a (c + 2ci)\v a \v a: (12a) 
T p = T + 2 Cl k a \v a \ 2 cos(a). (12b) 



A particular bisymmetric body is the sphere whose aero- 
dynamic characteristics (zero lift and constant drag co- 
efficient) are obtained by setting c\ = in (11). Elliptic- 
shaped bodies are also symmetric but, by contrast with 
the sphere, they do generate lift in addition to drag. The 
process of approximating measured aerodynamic char- 
acteristics with functions given by (11) is illustrated by 
the Figure 3 where we have used experimental data bor- 
rowed from [8, p. 19] for an elliptic-shaped body with 
Mach and Reynolds numbers equal to M = 6 and R e = 
7.96- 10 6 respectively. For this example, the identified co- 
efficients are Co = 0.43 and c\ = 0.462. Since missile-like 
devices are "almost" bisymmetric, approximating their 
aerodynamic coefficients with such functions can also be 
attempted. For instance, the approximation shown in 
Figure 4 has been obtained by using experimental data 
taken from [14, p. 54] for a missile moving at M = 0.7. 
In this case, the identified coefficients are Co = 0.1 and 
c\ = 11.55. In both cases, the match between experi- 
mental data and the approximating functions, although 
not perfect, should be sufficient for feedback control pur- 
poses. 

Other functions satisfying the condition (8) may also 
be of interest to model the aerodynamic characteristics 
of other, almost bisymmetric, bodies. For instance, by 
considering annular wings, functions of the form 

Cd{ol) = cq, Cx(a) = ci tan(a), 



The important result is the non-dependence of F p upon 
the angle of attack a, and thus on the vehicle's orienta- 
tion. The interest of this proposition is to point out the 
possibility of seeing a symmetric body subjected to both 
drag and lift forces as a sphere subjected to an equiva- 
lent drag force F p and powered by an equivalent thrust 
force T p = —T p k. The main condition is that the relation 
(8) must be satisfied. Obviously, this condition is com- 
patible with an infinite number of functions Cd and Cl- 
Let us point out a particular set of simple functions, al- 
ready considered in the 2D-case addressed in [13], which 
also satisfy the 7r-periodicity property w.r.t. the angle of 
attack a associated with bisymmetric bodies. 

Proposition 2 The functions Cd and Cl defined by 

(C D (a) = c + 2ci sin 2 (a) 
\ Cx(ce) = ci sin(2a) 



with Co and c\ denoting two real numbers, can be used to 
model classical linearly increasing lift and quasi-constant 
drag at small angles of attack, i.e. before the stall region 
is reached, in combination with functions (11) that are 
more representative of the physics for angles beyond the 
stall region. The possibility of forming in this way more 
complex functions that apply to variously shaped bodies 
and that (almost) satisfy the condition (8) over a large 
range of angles of attack is the subject of complementary 
studies. 



4 Control 

To illustrate the interest of the transformation evoked 
previously, let us consider the problem of stabilizing a 
reference velocity v r = (zb ? Jo? ko)x r = J, k)v r asymp- 
totically. The application of the control solution pro- 
posed in [6, Sec. III.D] to System (l)-(2), with (F a ,T) 



5 




two components of this vector are well defined. The in- 
terest of the invoked transformation, combined with (8), 
lies precisely there. As for the last component CJ3, since 
it does not influence the vehicle's longitudinal motion 
due to the symmetry about the axis Gk, it does not have 
to be defined at this point. This free degree of freedom 
can be used for complementary purposes involving, for 
instance, the angle /?. 

Let G (— 7r,7r] denote the angle between k and f p . In 
[6] , stability and convergence properties associated with 
the feedback control (13) are established by using the 
Lyapunov function candidate 



Figure 3. Aerodynamic coefficients of elliptic shaped bodies. 
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Figure 4. Aerodynamic coefficients of missile-like bodies. 

replaced by the equivalent drag force and thrust inten- 
sity (Fp,T p ) defined in Proposition 1, yields the follow- 
ing control expressions 



T = fa 3 + h\f p \v 3 , 



- k 2 \fp\v 2 
^2 = k 2 \f p \v 1 + 



f^S( ei )R T f p 



(\f P \ + U 2 \u? 

h\f P \f Pl f p r S(e 2 )R T f p 



with v := v 



(\f P \ + U 2 \u? 

v r , a r := j^v r (r , Jo, k )x r , 



(13a) 
, (13b) 

(13c) 



fa (?o, Jo,^o)/a (hJ,k)f a '= mg + F a - ma r , (14a) 
/p = (^o,Jo,^o)/p = (h3,k)f p := mg + F p - ma r , (14b) 

and fci,2,3 three positive real numbers. Note that, using 
(10b), the vector f p of coordinates of / p expressed in the 
fixed frame X is equal to mges — k a CD \v a \xa — rn%r, and 
is thus independent of the vehicle's orientation. There- 
fore, its time-derivative does not depend on the angular 
velocity vector uo and the above expressions of the first 



whose time-derivative along any trajectory of the con- 
trolled system is 



V 



-ki\f P \v 



k 3 , 2 

o — — tan 
k 2 



Assuming that v w and v r are bounded in norm up to 
their second time-derivatives, and provided that 3 S > 
such that \f p \ > 5, Vt G R + , one shows that the equilib- 
rium (v,0) = (0,0) of the controlled system is asymp- 
totically stable, with the domain of attraction equal to 

R 3 X (-7T,7T). 

In practice, the control law must be complemented 
with integral correction terms to compensate for al- 
most constant unmodeled additive perturbations. The 
solution proposed in [6] involves I v = (to,jo,ko)I v with 
I v := f Q x(s) ds, and x := Rv the longitudinal ve- 
locity error expressed in the inertial frame. Also, let 
h denote a smooth bounded strictly positive function 
defined on [0, +oo) satisfying the following properties 
([6, Sec. III.C]) for some positive constant numbers 77, /i, 
Vs G R, \h(s 2 )s\ < rj and < -§^(h(s 2 )s) < /jl. It then 
suffices to replace the definitions (14) of f a and f p by 



fa '= mg + F a -ma r + h(\I v \ 2 )I v 
fp := mg + F p - ma r + h(\I v \ 2 )I v 



(15a) 
(15b) 



in (13) to obtain a control which incorporates an inte- 
gral correction action and for which strong stability and 
convergence properties can also be proven (c.f. [6]). 

Remark 1 As in the case of velocity control, the 
position control presented in [6] can be adapted to 
stabilize reference trajectories. Videos of simulations 
with several reference trajectories can be found at 
|http://goo . gl/HKqtz| 
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5 Conclusion and perspectives 

The paper sets basic principles for the modeling and non- 
linear control of aerial vehicles subjected to strong aero- 
dynamic forces. Possible extensions are numerous. They 
concern in particular airplanes and other vehicles whose 
lift properties mostly rely on the use of large flat sur- 
faces (wings) which break the body symmetries here con- 
sidered. The dreaded stall phenomenon, when it is pro- 
nounced, is not compatible with transformations alike 
these evoked in the paper, because it forbids the exis- 
tence of equivalent drag forces that do not depend on the 
vehicle's angle of attack. Nor is it even compatible with 
the uniqueness of cruising equilibria and the objective of 
asymptotic stabilization, as pointed out in [12]. Its im- 
portance at both the modeling and control levels, and its 
consequences during transitions between hovering and 
lift-based-cruising, need to be studied and, if possible, 
attenuated via an adequate control design. Clearly, the 
control solution here proposed also calls for a multitude 
of complementary extensions and adaptations before it 
is implemented on a physical device. Let us just men- 
tion the production of the control torque allowing for 
desired angular velocity changes and the determination 
of corresponding low level control loops that take actu- 
ators' physical limitations into account -in relation, for 
instance, to the airspeed dependent control authority as- 
sociated with the use of flaps and rudders. Measurement 
and estimation of various physical variables involved in 
the calculation of the control law, other than the ever 
needed information about the vehicle's position and at- 
titude, such as the air velocity and the angle of attack, or 
the thrust force produced by a propeller, also involves a 
combination of hardware and software issues which are 
instrumental to implementation. 
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